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Abstract 

Background: MicroRNAs (miRNAs) are small, non-coding 21-25 nt RNA molecules that play an important role in 
regulating gene expression. Little is known about the expression profiles and functions of miRNAs in skin and their 
role in pigmentation. Alpacas have more than 22 natural coat colors, more than any other fiber producing species. 
To better understand the role of miRNAs in control of coat color we performed a comprehensive analysis of miRNA 
expression profiles in skin of white versus brown alpacas. 

Results: Two small RNA libraries from white alpaca (WA) and brown alpaca (BA) skin were sequenced with the aid 
of lllumina sequencing technology. 272 and 267 conserved miRNAs were obtained from the WA and BA skin 
libraries, respectively. Of these conserved miRNAs, 35 and 13 were more abundant in WA and BA skin, respectively. 
The targets of these miRNAs were predicted and grouped based on Gene Ontology and KEGG pathway analysis. 
Many predicted target genes for these miRNAs are involved in the melanogenesis pathway controlling 
pigmentation. In addition to the conserved miRNAs, we also obtained 22 potentially novel miRNAs from the WA 
and BA skin libraries. 

Conclusion: This study represents the first comprehensive survey of miRNAs expressed in skin of animals of 
different coat colors by deep sequencing analysis. We discovered a collection of miRNAs that are differentially 
expressed in WA and BA skin. The results suggest important potential functions of miRNAs in coat color regulation. 
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Background 

Since the first small RNA (sRNA) - Un4 was discovered 
in C. elegans, identification of small endogenous non- 
coding RNA (sncRNA) has received considerable atten- 
tion. SncRNAs include small interfering RNA [1], 
miRNA [2] and piwi-interacting RNA, which all regulate 
genes at the post- transcriptional level [3]. Among these 
sncRNAs, miRNAs are small, non-coding 21-25 nt RNA 
molecules that play an important role in regulating gene 
expression in animals and plants by promoting mRNA 
degradation and inhibiting mRNA translation [4]. How- 
ever, other studies revealed that some miRNAs may also 
function to induce gene expression [5]. Many miRNAs 
are evolutionary conserved in related species and some 
even show conservation between invertebrates and 
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vertebrates [6]. Generally, the functions of miRNAs are 
not limited to the regulation of developmentally timed 
events. Instead, they have diverse expression patterns 
and probably regulate many aspects of development and 
physiology [7]. To date, the miRBase database contains 
9072 mature miRNA products from 25 mammalian spe- 
cies. Recently, transcriptome profiling of miRNAs has 
led to identification of novel miRNAs, provided valuable 
information on tissue specific expression and evolution- 
ary conservation across species, and served as a corner- 
stone for subsequent functional analyses. 

In adult animals, both hair and skin color depend on 
pigment produced by melanocytes at the base of the epi- 
thelium [8]. Melanocytes in mammals and birds produce 
two types of melanin, black to brown eumelanin and yel- 
low to reddish brown pheomelanin [9,10]. The quality 
and ratio of eumelanin to pheomelanin influence the 
color of skin, fleece and eyes [9]. At present, a large 
number of genes have been found to affect hair and skin 
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color in humans and other vertebrate species [11-13], al- 
though the molecular and cellular mechanisms regulat- 
ing coat color in fiber-producing species, such as the 
alpaca, are not completely understood. 

Alpacas have a great variety of natural coat colors. 
Studies of the melanocortin-1 receptor (MC1R) and 
agouti signaling protein (ASIP) genes in alpacas of differ- 
ent coat colors revealed several polymorphisms that may 
be predictive of coat color [14-18]. Powell et al. reported 
that the M87V and S126G variants of MC1R may affect 
coat color with brown coat color animals showing 
homozygosity for the methionine and serine alleles [16]. 
Tyrosinase related protein 1 (TYRP1) is encoded by the 
brown locus and is involved in eumelanin synthesis [19]. 
Mutations in TYRP1 are associated with brown coat 
color in many species [20,21], but mutations in TYRP1 
associated with brown coat color in alpacas have not 
been reported to date. Previous studies in our laboratory 
revealed that the expression of many coat color genes 



including MC1R, TYRP1, tyrosinase (TYR), tyrosinase 
related protein 2 (TYRP2), nitric oxide synthase 2 (NOS2), 
paired box protein 3 (PAX3), and WNT3A was higher in 
skin of brown versus white alpacas, suggestive of a func- 
tional contribution to coat color regulation in this species. 
However, the potential contribution of miRNAs to dif- 
ferential expression (associated with coat color) of such 
genes has not been fully elucidated, primarily due to a 
lack of knowledge of the miRNA transcriptome profile 
in alpaca skin. 

Previous studies have addressed expression and function 
of miRNAs in mouse, goat and sheep skin [22-25]. For ex- 
ample, miR-203, a skin specific miRNA, was up-regulated 
in psoriatic plaques and shown to be involved in inflamma- 
tory responses and keratinocyte functions [26]. However, 
studies of the functional role of miRNAs in coat color regu- 
lation are limited. Dong et al. over-expressed miR-137 in 
mice and obtained a series of coat color changes in trans- 
genic animals [27]. Zhu et al. reported differential 
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Figure 1 Length distribution and abundance of unique sRNAs in white and brown alpaca skin. 
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expression of miR-25 in skin of white versus brown alpacas 
and a functional role of this miRNA in the regulation of ex- 
pression of microphthalmia-associated transcription factor 
(MITF), a transcription factor regulating expression of key 
genes involved in melanogenesis [28]. 

To better understand the potential role of miRNAs in 
the post-transcriptional regulation of genes linked to 
pigmentation, we have characterized the miRNA tran- 
scriptome in skin of alpaca of different coat colors by 
deep sequencing analysis. We report here the identifica- 
tion of hundreds of conserved miRNAs and candidate 
novel miRNAs in WA and BA skin. We also provide evi- 
dence of differential expression of miRNAs in WA and 
BA, foundational to future studies of the potential func- 
tional role for such miRNAs in post-transcriptional gene 
regulation linked to coat color. 

Results 

Overall complexity of sRNA pools between the WA and 
BA skin libraries 

To identify miRNAs expressed in alpaca skin, WA and 
BA skin small RNA libraries were analyzed by deep se- 
quencing. A total of 17,427,359 and 19,748,517 raw 
reads were obtained from WA and BA libraries respect- 
ively. After filtering out the low quality reads based on 
base quality value and removing adaptor sequences, 
16,834,090 (98.42%) and 19,069,600 (98.49%) clean reads 
were retrieved, corresponding to 1,058,335 and 977,308 
unique sequences, which were the associated counts of 
the identical sequence reads, respectively. The length 
distribution of sRNA reads is shown in Figure 1. The 
majority of the sRNAs were 21-24 nt in length, which is 
characteristic of sRNAs from metazoans. The most fre- 
quent length of sRNA sequenced was 22nt, which is 
identical to the classical size of Dicer cleavage products. 

The alpaca genome has been completely sequenced, 
but at low coverage and is not well annotated. Thus, the 
alpaca (http://www.ncbi.nlm.nih.gov) and Bos taurus 
genome databases (http://www.ncbi.nlm.nih.gov/nuccore/) 
were combined and used as a reference genome dataset 
for annotation of alpaca sRNAs. The total sRNA reads 
were mapped to the combined genome using the SOAP 
program, yielding 10,389,135 (61.71%) and 12,630,855 
(66.24%) genome matched reads in WA and BA librar- 
ies, respectively. In order to annotate the classes of all 
sRNAs, all clean reads were mapped using several data- 
bases, including GenBank, Rfam, animal miRNAs (miR- 
Base), exon and intron (Bos taurus) and repeat 
associated RNA. The representations of different types 
of sRNAs in WA and BA libraries are shown in Figure 2. 
By aligning the sRNAs to all known animal miRNA pre- 
cursors, 3820 (6,527,628 total reads) and 3823 
(7,583,101 total reads) unique sequences were identified 



as potentially conserved miRNA in WA and BA skin, 
respectively. 

Identification of conserved miRNAs in alpaca skin 

To further identify conserved miRNAs in alpacas, all 
sRNA sequences were mapped to known miRNAs in the 
miRBase 17.0 database. After alignment and additional 
sequence analysis, 272 and 267 individual candidate 
miRNAs from WA and BA skin libraries were identified, 
respectively. Approximately 155 miRNA families were 
represented in each library. Several miRNA families were 
represented by over 100,000 reads, such as Let-7, miR- 
146, miR-143, miR-21, miR-10, miR-378, miR-182, miR- 
8, and miR-30, which were the 9 most abundant miRNA 
families represented in WA and BA libraries (Figure 3). 
The read diversity of miRNA families indicated that dif- 
ferent miRNA families have different expression pat- 
terns. Besides the diversity in miRNA families, there was 
considerable variation in number of miRNA family 
members represented across families. For example, 13 
and 12 miR-154 family members were represented in 
WA and BA libraries, whereas many miRNA families 
only had 1 to 2 members (Figure 4). 

Differential expression of conserved miRNAs in white and 
brown alpaca skin 

In high throughput sequencing analysis, the number of 
reads is reflective of the expression level of miRNAs. 
Many conserved miRNAs exhibited dissimilar expression 
patterns in WA versus BA skin. There were 35 miRNAs 
showing markedly greater abundance in WA skin, and 
13 miRNAs more highly expressed in BA skin (Table 1). 
For example, miR-202, miR-542-5p, miR-424, miR-370 
and miR-22-3p were more abundantly expressed in WA 
skin, while miR-21 1, miR-184, miR-486, miR-885 and 
miR-451 were more highly expressed in BA skin. In 
order to validate the deep sequencing results, 10 miR- 
NAs were selected to determine their expression levels 
in white and brown alpaca skin by real time PCR. The 
results of real time PCR analysis confirmed the deep se- 
quencing results (Figure 5). For example, based on deep 
sequencing results, the expression level of miR-21 1 in 
WA was almost 20 times higher than observed in BA. 
Based on real time PCR results, the expression level of 
miR-21 1 in WA skin was also 13 times higher than 
observed in BA skin. 

Sequence divergence within known miRNAs 

The heterogeneity across sequenced reads at the 5' and 
3' ends from the same arms within individual known 
miRNAs detected were also analyzed. The 5' ends exhib- 
ited less heterogeneity than the 3' ends. For instance, 
miR-29c, miR-362-5p, and miR-222 were highly repre- 
sented and sequenced from the two libraries. The 
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predominant reads for these miRNAs obtained from 
deep sequencing data varied from the mature miRNA 
sequences in miRBase at the 3' end, and they should be 
regarded as the primary functional molecule (Figure 6). 
Collectively from the two libraries, the 5' and 3' end of 
109 alpaca miRNAs were different in sequence from cor- 
responding annotated miRNAs for each species in 



miRBase (Additional file 1). Furthermore, several abun- 
dant miRNAs with higher sequencing reads, like miR- 
30a, miR-324 and miR-330, were derived from different 
arms of the hairpin precursor than corresponding anno- 
tated miRNAs in miRBase (Figure 7). Thus, results sug- 
gest that the precursors of these miRNAs develop 
functional molecules from both arms in the alpaca. 
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Table 1 Significantly expressed miRIMAs in white and 



brown alpaca skin 


miR-name 


NEL in brown 


NEL in white 


FC 


White alpaca 


miR-202 


1 QQ77 
I .yyZ/ 


Qc: 993ZL 
yj. Z.Z.D'-t 


47 7861 


1 1 lln L \u-x 


U.I j/j 


I .OOJ J 


1 0 5741 
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0 3146 


2 8514 


9 0636 
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J. ZD 1 Z 


9c i Q7n 

ZD. 1 O/ U 


7 7ZL7P) 
/ ./H/ U 


miR-370 


0 5244 


J.OZ JO 


6 Q1 DO 
O.y I UU 


1 1 lln ZZ D|J 
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ZU/ .y/.Z.\j 


1 3^4. QAQQ 
I DD / -t.O\JyO 


O.J I oz 


1 1 lln Z/ a D\J 


71 6373 
/ I ,DjZj 


ZL9Q ZLR^A 


D.yyJ 1 


1 1 lln ^JU 


Z.KJy / D 


1 9 rmRQ 

I Z.UJO^ 


^ 7ZLRQ 


miR-^ZLr 
1 1 lln j^-L 


1 1 37QZL 


JO. 1 yJJ 


ZL Q3R3 


miR-1 11 
1 1 lln I ZZ 


^ Q7R1 
D.y / O I 


19. HQ77 
Ao.yJy / / 


ZL 7DD1 

H. /UU I 


miR-1 36 


U.J I H-O 


1 2475 


3 Q6S3 
J.^OJ J 


miR-4Q3 
1 1 im tyj 


0 891 5 


3 2672 


3 6648 


miR-?1 ^3 
1 1 lln Z 1 Od 


Z.^jDO 


1 D DQR6 
I u.u^oo 


3 ZL3RQ 


miR-671 


9 f)Q7fi 


0.07UO 


3 ?8M 
J.ZOJ I 


miR-4Q4 
1 1 lln tjt 


1 DZLRR 


J.UO^U 


1 QA^l 


miR-1 QDp 
1 1 lln 1 y\Ja 


4 2476 


1 1 7024 


7 7^1 

Z./ J J I 


miR-1 
1 1 im i od 


4 2476 


1 1 5836 


2 7271 


miR-1 94.3 
1 1 lln 1 Z^ta 


U./ ODD 


9 1 38^ 

Z. I JOJ 


9 71 87 
Z./ I o/ 


miR-1 Idh 
1 1 lln 1 zJ-x\J 


U./ ODD 


1 1 38^ 

Z. I JOJ 


2 71 87 


miR-1 4Ra 
1 1 II n 1 L \Oa 




A?6^ 

^l-ZOJ.^OOJ 


Z. JOU/ 


miR-?Dfi 
1 1 lln zuu 


Z./ ZDt? 


O.O^UO 


7 ^77D 

Z.JZ/ u 


miR-?(VL 
1 1 ii n zu^ 


1 8354 


46335 


2 5245 


miR-^AQ-^n 
1 1 lln jD7 op 


3 zl^i n 


R 7393 
O./jZj 


9 ^931 

Z.JZj I 


miR-^zm-^n 
1 1 lln D'HO op 


I .DZJD 


ZL D3QZL 


9 zLRzLQ 


miR-1 Q3a-3n 
1 1 lln 1 yJa D\J 


U.J / DO 


1 4257 


2 471 8 


miR-874 


0.8915 


2.1385 


2.3988 


miR-1 9a 


2.3073 


5.4057 


2.3429 


miR-1 43 


10818.3706 


24308.7093 


2.2461 


miR-1 296 


1 .4683 


3.2672 


2.2252 


miR-210 


3.6708 


7.9006 


2.1523 


miR-301a 


0.839 


1.7821 


2.1241 


miR-27a-5p 


55.009 


112.8662 


2.0518 


miR-1 46b 


25987.7501 


53201.0343 


2.0472 


miR-200a 


79.708 


162.1709 


2.0346 


miR-1 7-3p 


3.9854 


8.0788 


2.0271 


Brown alpaca 


miR-211 


47.4053 


2.9108 


16.2860 


miR-1 84 


232.0447 


22.2762 


10.4167 


miR-486 


983.9221 


151.6565 


6.4878 


miR-885 


5.1915 


0.9505 


5.4619 


miR-451 


1122.4147 


41 1 .0706 


2.7305 


miR-1 33a 


4.1427 


1.6039 


2.5829 


miR-299 


1 .7829 


0.7128 


2.5013 


miR-23a 


708.5623 


304.5606 


2.3265 


miR-25 


504.2581 


226.8017 


2.2233 



Table 1 Significantly expressed miRIMAs in white and 
brown alpaca skin (Continued) 



miR-23b-3p 


823.9292 


371.211 


2.2196 


miR-1 224 


6.6598 


3.0296 


2.1982 


miR-1 5b 


121.5547 


58.5716 


2.0753 


miR-1 93b 


207.1360 


100.4509 


2.0620 



Nel, normalized expression level of miRNA in a sample; FC, fold change of 
miRNAs between white and brown alpaca skin; Nel = actual miRNA count/total 
count of clean readsx 1000000. 



Target gene prediction for miRNAs in alpaca 

MiRNAs have functions in numerous developmental, 
physiological and pathological processes [29]. In order to 
provide insight into potential biological roles of miRNAs 
identified in alpaca skin, target prediction analysis was 
performed for alpaca skin miRNAs identified. The rules 
employed for target prediction were based on those 
reported by Allen et al. [30] and Schwab et al. [31]. 
There were 6317 genes predicted as the putative targets 
of conserved miRNAs in WA (272 miRNAs) and BA 
(267 miRNAs) skin, respectively. The predicted target 
genes were then subjected to gene ontology analysis for 
functional classification. Although the target functions 
of miRNA in WA and BA skin had some divergence, the 
top functions of targets were binding, transferase activ- 
ity, molecular transducer activity, and receptor activity 
(Additional file 2). The KEGG pathway was also used to 
further classify putative miRNA targets. All predicted 
targets were clustered into 257 pathways in WA and BA, 
respectively. Results were then interrogated for potential 
target genes implicated in coat color regulation and the 
melanogenesis pathway. Many genes in the melanogenesis 
pathway have already been shown to have a relationship with 
coat color in alpaca [32-35], such as pro-opiomelanocortin 
(POMC), melanocyte-stimulating hormone (MSH), MITF, 
kit oncogene (KIT) and TYR. Moreover, several miRNAs 
whose targets are involved in melanogenesis were 
expressed differentially in WA and BA skin. The target 
genes for the top 4 significantly expressed miRNAs in 
WA and BA included in the melanogenesis pathway are 
shown in Table 2. 



Identification of novel miRNAs in alpaca skin 

Because the alpaca genome has not been annotated, the 
bovine genome was utilized to identify potential novel 
miRNAs. It is well known that the important feature of 
miRNAs distinguishing them from other small RNAs is 
the ability of their genomic flanking sequences to fold 
into a hairpin structure [36]. The Mireap software was 
used to predict novel miRNAs using these properties. 
After filtering, 88 and 94 sequencing reads from WA 
and BA skin remained. Blast analysis against all mature 
miRNAs in miRBase, revealed 55 miRNAs from WA 
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Figure 5 Real time PCR validation of miRNA expression in skin of alpacas with white vs brown coat color. Abundance of miR-202, miR- 
21 1, miR-27a, miR-200a, miR-193b, miR-423, miR-484, miR-205, miR-137 and miR-145 were normalized relative to abundance of U6 snRNA. Bars in 
each panel represent the mean ± standard error (n=4 each), ** P < 0.05; *** P < 0.001. 



and 63 miRNAs from BA libraries that had no homology 
with other species. In these novel miRNAs, 31 miRNAs 
were common between WA and BA libraries. Thus, 
a total of 87 unique novel miRNAs were identified 
(Additional file 3). The hairpin structures for these novel 
miRNAs were denoted in Additional file 4. The lengths 
of these novel miRNAs ranged from 21 to 24nt, and the 
minimum free energy (MFE) of the novel miRNAs var- 
ied from -55.90 to -19.00 kcal mol 1 . The average MFE 
of novel miRNAs was -34.83 kcal mol" 1 , which is lower 



than tRNA and rRNA. This result was similar to 
results of Bonnet et al. [37] indicating that the major- 
ity of miRNAs exhibit a folding free energy which is 
lower than shuffled sequences. To further identify 
genuine miRNAs, the length of the RNAs was also 
considered. The MFE index (MFEI) was used to re- 
evaluate the potential novel miRNAs identified. Zhang 
et al. suggested that the MFEI for miRNAs was 0.97, 
significantly higher than tRNAs, rRNAs or mRNAs [38]. 
The MFEI of the 87 potential novel miRNAs ranged from 
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miR-29c 

ATCTCTTfiCACfiGGCTGACCGfiTTTCTCCTGGTGTTCAGAGTCTGTTTTTGTCTfiGCACCATTTGfiAftTCGGTTfiTGfiTGTAGGGGGA 

■(((((((((((...(((((((((((...(((((((((((( )))))...)))))))...))))))))))))).)))))))))• 



.CTAGCACCATTTGAAATCGGTT 


22 


2 


.CTAGCACCATTTGAAATCGGT 


21 


2 


. (TAGCAGCATTTGAAATCGG . [ 


19 


7(3 


. .TAGCACCATTTGAAATCGGTT 


21 


61 


. .TAGCACCATTTGAAATCGGTTA 


22 


49 


..TAGCACCATTTGAAATCGGT 


20 


42 


. .TAGCACCATTTGAAATCGGTTAT 


23 


3 


..TAGCACCATTTGAAATCG 


18 


1 


CACCATTTGAAATCGGTTA 


19 


1 



miR-362-5p 

CTCGAATCCTTGGAACCTAGGTGTGAGTGCTGTTCTAGTGCAACACAGCTATTCAAGGATTCAAA 



■ ■ (((((((((((■■■((((((((■(■■(( (■■■)))■■)■■)))))))))))))))))))■■■ length count 

. . |.AATCCTTGGAACCTAGGTGTGAGTGE| 26 249 

...AATCCTTGGAACCTAGGTGTGAG 23 50 

. . .AATCCTTGGAACCTAGGTGTGAGT 24 29 

. . .AATCCTTGGAACCTAGGTGTGA... 22 21 

. . .AATCCTTGGAACCTAGGTGTGAGTG 25 4 

. . .AATCCTTGGAACCTAGGTGTGAGTGCT 27 2 

ATCCTTGGAACCTAGGTGTGAG 22 1 

CTTGGAACCTAGGTGTGAGTGC 22 2 



miR-222 

GCTGCTGGAATGTGTAGGTACCCTCAATGGCTCAGTAGCCAGTGTAGATCCTGTCTTTTGTAATCAGTAGCTACATCTGGCTACTGGGTCTCTGATGGCATCTTCTACCT 

((...(...((((.((■(((..(((((((((((((((((((...((( )) )■■■)))))■))))))))))))))■■ )))■))■))))■)■■))■ length count 

I .AGCTACATCTCGCTACTGGGTcf] 23 316 

AGCTACATCTGGCTACTGGGTCTCT 25 219 

AGCTACATCTGGCTACTGGGTCTC 24 181 

.AGCTACATCTGGCTACTGGGT. 21 133 

AGCTACATCTGGCTACTGGG 20 74 

AGCTACATCTGGCTACTGGGTC 22 37 

AGCTACATCTGGCTACTGG 19 2 

AGCTACATCTGGCTACTGGGTCTCTGA 27 2 

CTACATCTGGCTACTGGGTCTCT 23 1 

CTACATCTGGCTACTGGGTCTC 22 1 

TACATCTGGCTACTGGGTCTC 21 1 

Figure 6 Heterogeneity at the 5' and 3' ends of the sequenced reads. Pre-miRNAs, structures, and multiple isoforms of expressed mature 
miR-29c, miR-362-5p, miR-222 sequences and their read counts are shown. The red underline sequences are annotated in miRBase17.0. The 
boxed sequences owned the highest read counts. 



0.38 to 2.45. Using the MFEI cutoff value of 0.97, 22 novel 
miRNAs present in the two libraries were identified 
(Table 3). Expression of these new miRNAs was typically 
several to a hundred counts lower than conserved miRNAs. 

Discussion 

Studies on expression and function of miRNAs in evolu- 
tionarily diverse agricultural species are limited. The al- 
paca is a domestic mammal specialized in fiber 
production [39] and has more than 22 natural coat col- 
ors [40-44]. In this report, we explored the miRNA tran- 
scriptome in alpaca skin using deep sequencing 
technology, and provided a foundation for future func- 
tional studies on the relationship between miRNA 



expression and development and function of skin includ- 
ing coat color. 

Many miRNAs have been identified by traditional San- 
ger sequencing and microarray methods. However, the 
high throughout-sequencing method which yields 
400,000 reads per run and identifies low abundance tags 
is a much more comprehensive method for miRNA pro- 
filing [45]. Using this sequencing approach, we found 22 
novel, as well as 272 and 267 conserved miRNAs in WA 
and BA skin libraries [2]. According to the new 
miRNA family classifications in Rfam (miRBase 17.0), 
272 miRNAs from WA and 267 miRNAs from BA 
were grouped into 155 families. Most of these identi- 
fied miRNA families are also conserved in other 
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niR-30a 

CTGTAAACATCCTCGACTGGAAGCTGTGAGGCTGCAGAAAGGCTTTCAGTCGGATGTTTGCAGCTGC 

((((((((((((..(((((((((((-(( ))-...))))))))))))))))))))))).... length count 



. TGTAAACATCCTCGACTGGAAG 




22 


1992 


.TGTAAACATCCTCGACTGGAAGC 




23 


425 


.TGTAAACATCCTCGACTGGAA 




21 


359 


.TGTAAACATCCTCGACTGG 




19 


189 


.TGTAAACATCCTCGACTGGA 




20 


119 


.TGTAAACATCCTCGACTGGAAGCT 




24 


12 




|.CTTTCAGTCGGATGTTTGCAGC J . . 


22 


11720 




CTTTCAGTCGGATGTTTGCAG 


21 


3718 




CTTTCAGTCGGATGTTTGCA 


20 


135 




CTTTCAGTCGGATGTTTGC 


19 


55 




CTTTCAGTCGGATGTTTGCAGCT. . 


23 


38 




TTTCAGTCGGATGTTTGCAGC . . . 


21 


5246 




TTTCAGTCGGATGTTTGCAG 


20 


1744 



miR-324 

AACTGGCTATGCCTCCCCGCATCCCCTAGGGCATTGGTGTAAAGCTGGAGACCCACTGCCCCAGGTGCTGCTGGGGGTTGTAGTCTGAC 

....((((((.-(((((.(((.(.(((.(((((.(((.a )))))-))))).))).).))).)))))..)))))).... len g th count 

CGCATCCCCTAGGGCATTGGTG 22 13 

CGCATCCCCTAGGGCATTGGT 21 7 

CGCATCCCCTAGGGCATTGGTGT 23 2 

22 35 



CCACTGCCCCAGGTGCTGCTGG 



.CCACTGCCCCAGGTGCTGCTG 21 7 

. .CACTGCCCCAGGTGCTGCTGG 21 1 

. . .ACTGCCCCAGGTGCTGCTGG 20 2 



miR-330 

CTTCGGCGATCACTGCCTCTCTGGGCCTGTGTCTTAGGCTCTGCAAGATCAACCGAGCAAAGCACACGGCCTGCAGAGAGGCAGCGCTCAGCTC 

(((((-(-(((((((((((((((((((-(((--(((( ))))-))))))).))))..))))))))))).).)).))). length count 

TCTCTGGGCCTGTGTCTTAGGCT 23 11 



.TCTCTGGGCCTGTGTCTTAG 20 9 

.TCTCTGGGCCTGTGTCTTAGG 21 6 

GCAAAGCACACGGCCTGCAGA 21 14 

GCAAAGCACACGGCCTGCAGAGA 23 9 

GCAAAGCACACGGCCTGCAGAG 22 8 

|.CAAAGCACACGGCCTGCAGAGA . [ 22 15 

CAAAGCACACGGCCTGCAGAGAGG 24 7 

CAAAGCACACGGCCTGCAGA 20 5 

CAAAGCACACGGCCTGCAGAG 21 5 

CAAAGCACACGGCCTGCAGAGAG 23 5 

AAAGCACACGGCCTGCAGAGA 21 7 



Figure 7 Annotated miRNA and maximal sequences from different arms of the precursors. The red underlined sequences are the 
annotated mature sequences in miRBase while the boxed sequences represent the maximum counts in sequencing database. 



mammalian species, such as mouse [23,46], human 
[26], goat and sheep [22]. Expression of these con- 
served alpaca miRNA families differed between WA 
and BA skin, based on both abundance and represen- 
tation within miRNA families. For example, miR-146 
and miR-8 families both had abundant reads in WA 



and BA skin libraries, but miR-146b (miR-146 family) 
and miR-141 (miR-8 family) were of greater abundance 
in WA than BA. This significant variation maybe attribu- 
ted to the timing and location of miRNA expression dur- 
ing skin development in adult alpacas. In many miRNA 
families, there was a predominant member whose reads 



Table 2 Targets for top 4 significantly expressed miRNAs in white and brown alpaca involved in melanogenesis 

Name Predicted targets involved in melanogenesis 

miR-21 1 CPNE8, PRKCG, TMEM135, ARMC8, DCT, PLCB1,WNT3 

miR-424* GNAI2, TCF7, CREB3, FZD7, MYOF, PRKACB, KIT, GNAI1, PLCB1, SYT1 1, SYT1, MCTP1, FZD1, SYT4, CALM 

miR-202 CPNE8,TCF7, KRAS, WNT7A, PRKCB, KIT, CAMK2D, CABP2, EDNRB, WNT5A, FZD4 

miR-184 CAMK2A, FAM193B, MYOF, ADCY5, ADCY8, PRKCB, PLCB3, RAF1, SYT1, WNT1 1, PRKACA, WNT16, PLCB2 



Tian et al. BMC Genomics 2012, 13:555 
http://www.biomedcentral.eom/1 471 -21 64/1 3/555 



Page 9 of 12 



Table 3 Novel miRNAs identified in alpaca skin 



Name 


Sequence 


MFEI 


Alpaca-novel-82 


TCGCTCTCCTGCTCGCTCTGC 


0.9789 


Alpaca-novel-60 


TCTGGAGGACGCCGCTCGCGCTC 


0.9811 


Alpaca-novel-56 


GAATGGTGCTCCCTGGAATOT 


0.9974 


Alpaca-novel-24 


AGGGAACC1TGAAAAGCTGAAG 


0.9986 


Alpaca-novel-83 


CAGGACCTGGGGACACCATOT 


1 .0079 


Alpaca-novel-37 


TGmCTCAGAAGACTGTAGT 


1.0155 


Alpaca-novel-17 


TGAACGGGGCCCTOTGGTAG 


1 .0246 


Alpaca-novel-61 


CCTCACCTGTCATOTCCCAGA 


1 .0324 


Alpaca-novel-31 


ATOGCATGTCCTGGAATGAG 


1 .0554 


Alpaca-novel-53 


CCAAACCAGTOTGCCTGTAG 


1 .0602 


Alpaca-novel-45 


ACmGGAmGAGTCTCTGGT 


1.1026 


Alpaca-novel-8 


ACAAAGTOGCTGCCTGTGAGC 


1.1197 


Alpaca-novel-9 


TOGTOAATACGTGGATGTGG 


1.1877 


Alpaca-novel-33 


ATCCGTGGTOGTOAATATGC 


1 .2365 


Alpaca-novel-32 


TGAATGGCACCTOTGAGTAGA 


1.2418 


Alpaca-novel-76 


AACTCAGTGTCAGATAGGAAGA 


1 .2434 


Alpaca-novel-27 


TOGTGATGTOGTOAAAGAG 


1.3272 


Alpaca-novel-44 


TCTGAAAmAAATGTAACCGG 


1 .3920 


Alpaca-novel-43 


TOTOTOGTOAATAGTA^ 


1 .5895 


Alpaca-novel-68 


TAATAGGTGATCAAATGAATGA 


1.9358 


Alpaca-novel-81 


ATOATCmGACTATAACTG 


2.0412 


Alpaca-novel-8 


ACAAAGTOGCTGCCTGTGAGC 


2.4490 



MFEI, MFE index. 



were significantly higher than other members. For ex- 
ample, the Let- 7 miRNA family contains 10 members. 
The reads for let-7a were more than one million, but let- 
7i only had five thousand reads. It is possible that the 
predominant member in each miRNA was responsible 
for the regulatory role at the developmental stage in 
which skin samples were collected in the present studies. 

We compared the expression of conserved miRNAs 
between WA and BA skin and identified differentially 
expressed miRNAs from the two libraries. The expres- 
sion of miRNAs has an established relationship with cell 
proliferation and differentiation. As such, differences in 
miRNA expression between animals possibly reflect 
alterations in growth and differentiation. Heterochrony, 
or change in the timing of developmental growth, is 
known to alter the overall shape of developing tissues, 
leading to differences in physiology [47]. Alternatively, 
there are pronounced phenotypic differences between 
WA and BA which include fiber diameter, staple length, 
hair color and other hair characteristics [48,49]. Thus, 
differences in miRNA abundance may also be related to 
hair color and other fiber characteristics. The functional 
role of such differentially expressed miRNAs is the focus 
of future investigation. 

In animals, miRNAs target mRNAs with partial se- 
quence complementarity [50]. One miRNA is inferred to 



influence expression of hundreds of mRNAs via this 
flexible recognition. There were many differentially 
expressed miRNAs identified between WA and BA skin, 
so the next step was to search for putative specific genes 
targeted by these miRNAs. The major physiological dif- 
ference of interest between the samples analyzed was 
coat color, so we focused on potential target genes regu- 
lating pigmentation. Based on KEGG pathway analysis, 
we found 117 potential miRNA target genes involved in 
melanogenesis. Many miRNAs targeting these genes 
have been reported. For examples, miR-137 [51] and 
miR-25 [28] were previously shown to regulate the ex- 
pression of MITF, which is an important regulator of 
melanocyte growth, maturation, apoptosis and melano- 
genesis. In the present study, miR-137 and miR-25 were 
found in both libraries and their target genes also 
included MITF. Haflidadottir et al. showed that miR-148 
affects Mitf mRNA expression in melanoma cells 
through a conserved binding site in the S'-UTR 
sequences of mouse and human Mitf [52]. In our stud- 
ies, miR-148 showed greater expression in WA versus 
BA. So, it is possible that miR-148 may contribute to 
coat color differences between WA and BA. 

The present studies also revealed that both miR-193a 
and miR-193b are differentially expressed in WA and BA 
skin. Gao et al. reported that Kit expression was subject to 
post- transcriptional regulation by miR-193a and miR-193b 
in human acute myeloid leukemia [53,54]. In humans, Kit 
contributes to the regulation of skin pigmentation [55] 
and we have previously reported Kit expression is higher 
in skin of BA versus WA [34]. Thus, miR-193a and miR- 
193b may play a role in coat color determination in alpa- 
cas via regulating the expression of Kit. 

Conclusions 

Studies reported provide novel information on miRNA 
expression in alpaca species, including conserved miR- 
NAs and potential novel miRNAs, and insight into dif- 
ferential miRNA expression associated with coat color 
and respective target genes. The predicted targets of the 
conserved miRNAs include 117 genes related to melano- 
genesis. Further validation of the functional role of dif- 
ferentially expressed miRNAs and their targets will 
provide better understanding of the miRNA regulatory 
mechanisms influencing coat color. 

Methods 

Sample collection and total RNA extraction 

Housing and care of alpacas and collection of skin sam- 
ples for use in described experiments were conducted in 
accordance with the International Guiding Principles for 
Biomedical Research Involving Animals (http://www. 
cioms.ch/frame 1985 texts of guidelines, html). Three 
white and three brown healthy adult alpaca (lama paco) 
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were selected from the alpaca farm in Shanxi Agricul- 
ture University (Shanxi, China) for sample collection. 
Three pieces of skin (2cm x3 cm) from the hindquarter 
were collected under local anesthesia and immediately 
stored in liquid nitrogen. Total RNA from the samples 
was extracted using Trizol reagent (Invitrogen, Carlsbad, 
CA, USA) according to the manufacturer s instructions. 
The RNA integrity was evaluated by gel electrophoresis 
and the RNA concentration was measured by OD260 
reading using a Nanodrop spectrophotometer. 

Small RNA library preparation and sequencing 

For deep sequencing, total RNA samples were size- 
fractionated on a 15% PAGE gel, and the 16-30 nt frac- 
tion was collected. The 5'RNA adaptor and 3'RNA 
adaptor were ligated to the RNA pools, and the RNAs of 
64-99 nts were isolated through gel elution and ethanol 
precipitation. PCR products were purified and small 
RNA libraries were sequenced using Illumina Genome 
Analyzer. Sequencing was carried out at the ShenZheng 
Genomics Institute (BGI, China). 

Initial processing of reads and identification of conserved 
miRNAs 

Raw sequence reads were produced by the Illumina 
HiSeq2000 and processed into clean full length reads by 
the BGI small RNA pipeline. During the procedure, all low 
quality reads, including 3' and 5' adapter contaminants 
were removed. The remaining high quality sequences lar- 
ger than 30nt and smaller than 18nt were discarded. Then, 
all high quality sequences were mapped to Vicugna pacos 
genome (http://www.ncbi.nlm.nih.gov/nuccore/206581 138) 
and Bos Taurus genome (http://hgdownload.cse.ucsc.edu/ 
goldenPath/bosTau4/bigZips/bosTau4.fa.gz) by SOAP [56]. 
All clean reads were compared with the sequences of 
non-coding RNA (rRNA, tRNA, snRNA, snoRNA) avail- 
able in Rfam (http://www.sanger.ac.uk/software/Rfam) 
[57] and the GenBank non-coding RNA database (http:// 
www.ncbi.nlm.nih.gov/) to classify degradation fragments 
of non-coding RNA. We followed the rule [58]: rRNAetc, 
in which Genbank > Rfam> known miRNA > repeat > 
exon > intron. Subsequently, the remaining sequences 
were aligned to the whole animal miRNAs deposited in 
miRBase 17.0 to obtain the known miRNAs in alpaca as 
well as base bias at the first position of identified miRNAs 
and at each position of all identified miRNAs respectively. 
If a known miRNA in alpaca showed homology with less 
than two mismatches (or >90% homology), it was consid- 
ered as evolutionary conserved. 

Validation of conserved miRNA expression 

Ten conserved miRNAs were selected for validation by 
quantitative real time PCR analysis. Total RNA was iso- 
lated from eight alpaca skin samples (4 brown, 4 white). 



One ug of DNase-treated RNA was converted to cDNA 
using specific stem-loop RT primer (Additional file 5) 
and HiFi-MMLV cDNA kit mix (CWBIO, Beijing, 
China). The cDNA was then used for real time PCR 
quantification of miRNAs using miRNA specific primers 
(Additional file 5) in combination with the universal pri- 
mer. U6 snRNA was used as an endogenous control 
(primers listed in Additional file 5). Quantitative real- 
time PCR was performed in triplicate on the Stratagene 
Mx3005P system. The 25 uL PCR reaction included 
12.5uL SYBR Premix Ex Taq™ II (TaKaRa, Dalian, 
China), 1.0 uL specific forward primer, luL universal primer, 
0.5 uL ROX reference dye, 2.5 uL diluted (4 times) cDNA 
and 7.5 uL water. Cycling parameters were 95°C for 30 sec 
followed by 40 cycles of 95°C for 5 sec, 56 °C or 58 °C 
for 20 sec and 72 °C for 15 sec. Melting curve analyses 
were performed following amplifications. Quantification 
of selected miRNA transcript abundance was performed 
using the comparative threshold cycle (CT) method 
established by Livak et al. [59]. The abundance of 
selected miRNAs was normalized relative to that of U6 
snRNA [60]. 

Identification of novel miRNAs 

The prediction of alpaca candidate novel miRNAs was 
conducted using criteria that were previously developed 
for animal miRNA prediction [61]. The criteria were as 
follows: the putative mature sequence must reside at the 
stem region and its size was limited to 20-24nt; the fre- 
quency of putative mature sequence should not be below 
5; the folding free energy of the stem-loop structure was 
limited to below -18 kcal mol 1 ; the maximum tolerance 
of a bugle size was 4 nucleotides; the maximum size of 
difference between miRNA and miRNA* was 5 nucleo- 
tides and 35 nucleotides; the sequence asymmetry be- 
tween miRNA and miRNA* could not exceed 5 
nucleotides. The selected sequences were then folded 
into a secondary structure using mFold3.2 (http://mfold. 
rit.albany.edu). If a perfect stem-loop structure was 
formed, the small RNA sequence was located at one arm 
of the stem and met the above criteria, the small RNA 
was then considered as a candidate novel alpaca miRNA. 
The formula for calculating MFEI was: MFEI = [(MFE / 
length of the RNA sequence) xlOO] / (G + C )%]. 

Target gene prediction 

Mireap software (https://sourceforge.net/projects/mireap/) 
was used to predict the target genes for conserved miRNAs 
identified. The criteria were as follows: no more than four 
mismatches between miRNA & target (G-U bases count as 
0.5 mismatches); no more than two adjacent mismatches 
in positions 2-12 of the miRNA/target duplex (5' of 
miRNA); no mismatches in positions 10-11 of miRNA/ 
target duplex; no more than 2.5 mismatches in positions 
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1-12 of the miRNA/target duplex (5' of miRNA); MFE 
of the miRNA/target duplex should be >75% of the 
miRNA bound to its perfect complement. 

Additional files 



Abbreviations 

miRNAs: MicroRNAs; WA: White alpaca; BA: Brown alpaca; sRNA: Small RNA; 
sncRNA: Small endogenous non-coding RNA; MC1R: Melanocortin-1 receptor; 
ASIP: Agouti signaling protein; TYRP1: Tyrosinase related protein 1; 
TYR: Tyrosinase; T/RP2: Tyrosinase related protein 2; N0S2: Nitric oxide 
synthase 2; PAX3: Paired box protein 3; MFE: Minimum free energy; 
MFEI: MFE index; MITF: Microphthalmia-associated transcription factor; 
POMC: Pro-opiomelanocortin; MSH: Melanocyte-stimulating hormones; 
KIT: Kit oncogene. 

Competing interests 

The authors declare that they have no competing interests. 
Authors' contributions 

XT designed the study, performed the analysis, analyzed the data and wrote 
the paper. J BY and GWS participated in the design of the study, analyzed the 
data and critically revised the manuscript. DLZ and YFS took part in 
collecting samples and the data analysis. JBJ, RWF, HDW, XLM, XYH, JPH, 
HQL, JJG and XJY conceived of the study and critically revised the 
manuscript. CSD participated in the design of the study, analyzed the data 
and participated in writing. All authors read and approved the final 
manuscript. 

Acknowledgements 

We thank Dr. Xinmin Li (Department of Pathology & Laboratory Medicine, 
David Geffer School of Medicine, UCLA), professor Lihua Lu, Yun Zhang for 
useful advice, discussion, and assistance and Mr. Wang Ying (China alpaca 
farm, Zhuangzi, Yuci, China) for providing experimental animals. The study 
was sponsored by the National Natural Science Foundation of China 
(30571070, 2006031060). 

Author details 

College of Animal Science and Technology, Shanxi Agricultural University, 
Taigu 030801, People's Republic of China. 2 Division of Animal and Nutritional 
Sciences, Laboratory of Animal Biotechnology and Genomics, West Virginia 
University, Morgantown, WV 26506, USA. laboratory of Mammalian 
Reproductive Biology and Genomics, Departments of Animal Science and 
Physiology, Michigan State University, East Lansing, Ml 48824, USA. 

Received: 9 April 2012 Accepted: 8 October 2012 
Published: 16 October 2012 



References 

1. Fire A, Xu SQ, Montgomery MK, Kostas SA, Driver SE, Mello CC: Potent and 
specific genetic interference by double-stranded RNA in Caenorhabditis 
elegans. Nature 1998, 391 (6669) :806-81 1. 

2. Carthew RW, Sontheimer EJ: Origins and Mechanisms of miRNAs and 
siRNAs. Cell 2009, 136(4):642-655. 

3. Grivna ST, Beyret E, Wang Z, Lin H: A novel class of small RNAs in mouse 
spermatogenic cells. Gene Dev 2006, 20(1 3):1 709-1 714. 

4. Bartel DP: MicroRNAs: Genomics, biogenesis, mechanism, and function. 
Ce//2004, 116(2):28 1-297. 

5. Li LC, Pookot D, Noonan EJ, Dahiya R: MicroRNA-373 induces expression 
of genes with complementary promoter sequences. P Natl Acad Sci Usa 
2008, 105(5):1608. 

6. Landgraf P, Rusu M, Sheridan R, Sewer A, lovino N, Aravin A, Pfeffer S, Rice 
A, Kamphorst AO, Landthaler M: A Mammalian microRNA expression atlas 
based on small RNA library sequencing. Cell 2007, 1 29(7):1401 -1414. 

7. He L, Hannon GJ: MicroRNAs: Small RNAs with a big role in gene 
regulation. Nat Rev Genet 2004, 5(7)522-531. 

8. Barsh G, Cotsarelis G: How hair gets its pigment. Cell 2007, 1 30(5):779-781 . 

9. Ito S, Wakamatsu K, Ozeki H: Chemical analysis of melanins and its 
application to the study of the regulation of melanogenesis. Pigm Cell Res 
2000, 13:103-109. 

10. Ito S, Wakamatsu K: Chemistry of mixed melanogenesis — pivotal roles of 
dopaquinone. Photochem Photobiol 2008, 84(3):582-592. 

11. Lamason RL: SLC24A5, a putative cation exchanger, affects pigmentation 
in zebrafish and humans. Science 2005, 31 0(5755):1 782-1 786. 

12. Passeron T, Valencia JC, Bertolotto C, Hoashi T, Le Pape E, Takahashi K, 
Ballotti R, Hearing VJ: SOX9 is a key player in ultraviolet B-induced 
melanocyte differentiation and pigmentation. P Natl Acad Sci Usa 2007, 
104(35):1 3984-1 3989. 

13. Chintala S: Slc7a11 gene controls production of pheomelanin pigment 
and proliferation of cultured cells. P Natl Acad Sci Usa 2005, 102 
(31):10964-10969. 

14. Feeley N, Munyard K: Characterisation of the melanocortin-1 receptor 
gene in alpaca and identification of possible markers associated with 
phenotypic variations in colour. Anim Prod Sci 2009, 49(8):675-681. 

15. Bathrachalam C, Manna VL, Renieri C, Terza AL: Asip and MCIR cDNA 
polymorphism in alpaca. In Fibre production in South American camelids 
and other fibre animals; 2012:93-96. 

16. Powell AJ, Moss MJ, Tree LT, Roeder BL, Carleton CL, Campbell E, Kooyman 
DL: Characterization of the effect of melanocortin 1 receptor, a member 
of the hair color genetic locus, in alpaca (Lama pacos) fleece color 
differentiation. Small Ruminant Res 2008, 79(2-3):1 83-187. 

1 7. Hart KW, Barker SJ, Skerritt JW, Vercoe PE: The gene for white fleece in 
alpacas is homologous to AWt, the gene for white fleece in sheep. In 50 
years of DNA: Proceedings of the Fifteenth Conference, Association for the 
Advancement of Animal Breeding and Genetics, Melbourne, Australia, 7-11 July 
2003. 2003:82-85. 

18. Feeley N, Bottomley S, Munyard K: Three novel mutations in ASIP 
associated with black. J Agr Sci 201 1, 149:529-538. 

19. Zdarsky E, Favor J, Jackson IJ: The molecular basis of brown, an old mouse 
mutation, and of an induced revertant to wild type. Genetics 1990, 
126(2):443-449. 

20. Rieder S, Taourit S, Mariat D, Langlois B, Guerin G: Mutations in the agouti 
(ASIP), the extension (MC1R), and the brown (TYRP1) loci and their 
association to coat color phenotypes in horses (Equus caballus). 
Mamm Genome 2001, 12(6)450-455. 

21. Schmidt-Kuntzel A, Eizirik E, O'Brien SJ, Menotti-Raymond M: Tyrosinase 
and tyrosinase related protein 1 alleles specify domestic cat coat 
color phenotypes of the albino and brown loci. J Hered 2005, 
96(4):289-301. 

22. Zhang WG, Wu JQ, Li JQ, Midori Y: A subset of skin-expressed microRNAs 
with possible roles in goat and sheep hair growth based on expression 
profiling of mammalian microRNAs. OMICS: A Journal of Integrative Biology 
2007, 11(4):385-396. 

23. Andl T, Murchison EP, Liu F, Zhang Y, Yunta-Gonzalez M, Tobias JW, Andl 
CD, Seykora JT, Hannon GJ, Millar SE: The miRNA-processing enzyme dicer 
Is essential for the morphogenesis and maintenance of hair follicles. 
Curr Biol 2006, 1 6(1 0)4041-1049. 

24. Bostjancic E, Glavac D: Importance of microRNAs in skin morphogenesis 
and diseases. Acta Dermatovenerol Alp Panonica Adriat 2008, 17(3):95. 



Additional file 1: Table SI. Heterogeneity at the 5' and 3' ends of the 
sequenced reads. The list of alpaca miRNAs different with miRNAs in 
miRBase based on sequencing counts. 

Additional file 2: Table S2. Gene Ontology analysis of miRNAs targets. 
Targets of conserved miRNAs were grouped by Gene Ontology. 

Additional file 3: Table S3. Novel miRNAs in white and brown alpaca 
skin. Sequences and MFE of candidate novel miRNAs in white and brown 
alpaca. 

Additional file 4: Structure S4. Stem-loop structure of novel miRNAs. 
The stem-loop structure of 87 novel miRNAs in white alpaca and brown 
alpaca. The red underlined sequences are the predicted mature 
sequences. 

Additional file 5: Table S5. miRNAs primers for real time PCR. List of 
miR-202, miR-211, miR-27a, miR-200a, miR4 93b, miR-423 miR-484, miR- 
205, miR-137, miR-145 and U6 specific stem-loop RT primer, specific 
primers and universal primer. 



Tian et al. BMC Genomics 201 2, 1 3:555 Page 1 2 of 1 2 

http://www.biomedcentral.eom/1 471 -21 64/1 3/555 



25. Yi R, Poy MN, Stoffel M, Fuchs E: A skin microRNA promotes 
differentiation by repressing 'sternness'. Nature 2008, 452(71 84):225-229. 

26. Sonkoly E, Wei T, Janson PCJ, Saaf A, Lundeberg L, Tengvall-Linder M, 
Norstedt G, Alenius H, Homey B, Scheynius A: MicroRNAs: novel regulators 
involved in the pathogenesis of psoriasis? PLoS One 2007, 2(7):e610. 

27. Dong C, Wang H, Xue L, Dong Y, Yang L, Fan R, Yu X, Tian X, Ma S, Smith 
GW: Coat color determination by miR-137 mediated down-regulation of 
microphthalmia-associated transcription factor in a mouse model. 
RNA 2012, 18:1679-1686. 

28. Zhu Z, He J, Jia X, Jiang J, Bai R, Yu X, Lv L, Fan R, He X, Geng J: MicroRNA- 
25 functions in regulation of pigmentation by targeting the transcription 
factor MITF in alpaca (Lama pacos) skin melanocytes. Domest Anim 
Endocrin 2010, 38(3)200-209. 

29. Sand M, Gambichler T, Sand D, Skrygan M, Altmeyer P, Bechara F: 
MicroRNAs and the skin: Tiny players in the body's largest organ. J Dermatol 
Sci 2009, 53(3):169-175. 

30. Allen E, Xie Z, Gustafson AM, Carrington JC: microRNA-directed phasing 
during trans-acting siRNA biogenesis in plants. Cell 2005, 121(2):207-221. 

31. Schwab R, Palatnik JF, Riester M, Schommer C, Schmid M, Weigel D: 
Specific effects of microRNAs on the plant transcriptome. Dev Cell 2005, 
8(4):5 17-527. 

32. Jiang J, Yu X, He J, Dong C: Study on the expression and localization of 
tyrosinase in skin of alpacas with different coat colors. Scientia Agricultura 
Scinica 2010, 43(12):2555-2560. 

33. Dong Y, Cao J, Wang H, Zhang J, Zhu Z, Bai R, Hao HQ, He X, Fan R, Dong 
C: Nitric oxide enhances the sensitivity of alpaca melanocytes to 
respond to a-melanocyte-stimulating hormone by up-regulating 
melanocortin-1 receptor. Biochem Bioph Res Co 2010, 396(4):849-853. 

34. Jiang J, Yu XJ, Tian X, He JP, Dong CS: SCF/c-KIT signal regulates the 
proliferation and division of alpaca (lama pacos) hair follicle melanocytes 
and its mechanism. Acta Anatomica Sinica 201 1, 42(1 ):99-1 03. 

35. Ren Y, Ren B, Fan R, Zhu Z, Yang Y, Li H, Dong C: Cloning and Expression 
Level Analysis of MC1R Gene in Alpacas with Different Coat Color. 
Acta Veterinaria et Zootechnica Sinica 2009, 40(6):5. 

36. Lagos-Quintana M, Rauhut R, Lendeckel W, Tuschl T: Identification of novel 
genes coding for small expressed RNAs. Science 2001, 294(5543):853. 

37. Bonnet E, Wuyts J, Rouze P, Van de Peer Y: Evidence that microRNA 
precursors, unlike other non-coding RNAs, have lower folding free 
energies than random sequences. Bioinformatics 2004, 20(1 7):291 1-2917. 

38. Zhang B, Pan X, Cox S, Cobb G, Anderson T: Evidence that miRNAs are 
different from other RNAs. Cell Mol Life Sci 2006, 63(2):246-254. 

39. Cecchi T, Valbonesi A, Passamonti P, Gonzales M, Antonini M, Renieri C: 
Quantitative variation of melanins in alpaca (Lama pacos L.). ltd J Anim 
Sci 2011, 10(3):e30. 

40. McGregor B: Production, attributes and relative value of alpaca fleeces in 
southern Australia and implications for industry development. Small Ruminant 
/?es 2006, 61 (2-3):93-1 1 1 . 

41. Hoffman E: The complete alpaca book, 2nd edn. Santa Cruz, California: Bonny 
Doon Press; 2006. 

42. Lauvergne JJ, Marinez G, Ayala C, Rodriguez T: Supreme-Project: 
Identification of a primary population of South American domestic 
camelids in the province of Antonio Quijarro and Enrique Baldivieso 
(Department of Botosi, Bolivia) using the phenotypic variation of coat 
colour. In Progress in South American camelids Reasearch. Edited by Gerken 
M, Renieri C. The Netherlands: Wageningen Pers; 2001:64-74. 

43. Bartolome E, Sanchez MJ, Peha F, Horcada A: Testing objective metrics for 
the differentiation of coat colours in a Spanish alpaca population. In 
Fibre production in South American camelids and other fibre animals. Edited 
by Perez-Cabal MA, Gutierrez J, Cervantes I, Alcalde MJ. The Netherlands: 
Wageningen; 2011:79-86. 

44. Fan R, Yang G, Dong C: Study of hair melanin in various hair color Alpaca 
(Lama pacos). Asian Austral J Anim Sci 2010, 23:444-449. 

45. Shendure J, Ji H: Next-generation DNA sequencing. Nat Biotechnol 2008, 
26(1 0):1 135-1 145. 

46. Yi R, O'Carroll D, Pasolli HA, Zhang Z, Dietrich FS, Tarakhovsky A, Fuchs E: 
Morphogensis in skin is governed by a discrete sets of differentially 
expressed microRNAs. Nat Genet 2006, 38(3):356-362. 

47. Ason B, Darnell DK, Wittbrodt B, Berezikov E, Kloosterman WP, Wittbrodt 
J, Antin PB, Plasterk RHA: Differences in vertebrate microRNA 
expression. P Natl Acad Sci Usa 2006, 103(39):1 4385-1 4389. 



48. McGregor BA, Butler KL: Sources of variation in fibre diameter attributes 
of Australian alpacas and implications for fleece evaluation and animal 
selection. Aust J Agr Res 2004, 55:433-442. 

49. Lupton C, McColl A, Stobart R: Fiber characteristics of the Huacaya alpaca. 

Small Ruminant Res 2006, 64(3):21 1-224. 

50. Karginov FV, Conaco C, Xuan Z, Schmidt BH, Parker JS, Mandel G, Hannon 
GJ: A biochemical approach to identifying microRNA targets. P Natl Acad 
Sci Usa 2007, 1 04(49): 1 9291 -19296. 

51. Bemis LT, Chen R, Amato CM, Classen EH, Robinson SE, Coffey DG, Erickson 
PF, Shellman YG, Robinson WA: MicroRNA- 137 targets microphthalmia- 
associated transcription factor in melanoma cell lines. Cancer Res 2008, 
68(5):1 362-1 368. 

52. Haflidadottir BS, Bergsteinsdottir K, Praetorius C, Steingnmsson E: miR-148 
regulates Mitf in melanoma cells. PLoS One 2010, 5(7):e1 1574. 

53. Gao X, Lin J, Gao L, Li Y, Wang L, Yu L: MicroRNA-193b regulates c-Kit 
proto-oncogene and represses cell proliferation in acute myeloid 
leukemia. Leukemia Res 201 1, 35(9):1 226-1 232. 

54. Gao X, Lin J, Li Y, Gao L, Wang X, Wang W, Kang H, Yan G, Wang L, Yu L: 
MicroRNA-193a represses c-kit expression and functions as a 
methylation-silenced tumor suppressor in acute myeloid leukemia. 
Oncogene 201 1, 30(3 1):341 6-3428. 

55. Kasamatsu S, Hachiya A, Higuchi K, Ohuchi A, Kitahara T, Boissy RE: 
Production of the soluble form of KIT, s-KIT, abolishes stem cell factor- 
induced melanogenesis in human melanocytes. J Invest Dermatol 2008, 
128(7):1 763-1 772. 

56. Li R, Li Y, Kristiansen K, Wang J: SOAP: short oligonucleotide alignment 
program. Bioinformatics 2008, 24(5):71 3-714. 

57. Griffiths-Jones S, Moxon S, Marshall M, Khanna A, Eddy SR, Bateman A: 
Rfam: annotating non-coding RNAs in complete genomes. Nucleic Acids 
Res 2005, 33(suppl 1 ):D1 21 -D1 24. 

58. Calabrese JM, Seila AC, Yeo GW, Sharp PA: RNA sequence analysis defines 
Dicer's role in mouse embryonic stem cells. P Natl Acad Sci Usa 2007, 
104(46):1 8097-1 81 02. 

59. Livak KJ, Schmittgen TD: Analysis of relative gene expression data using 
real-time quantitative PCR and the 2~ AACT method. Methods 2001, 
25(4):402-408. 

60. Choong ML, Yang HH, McNiece I: MicroRNA expression profiling during 
human cord blood-derived CD34 cell erythropoiesis. Exp Hematol 2007, 
35(4):55 1-564. 

61. Berezikov E, Cuppen E, Plasterk RHA: Approaches to microRNA discovery. 

Nat Genet 2006, 38(6s):S2-S7. 



doi:1 0.1 1 86/1 471 -21 64-1 3-555 

Cite this article as: Tian et al.: Identification and characterization of 
microRNAs in white and brown alpaca skin. BMC Genomics 2012 13:555. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at 
www.biomedcentral.com/submit 



o 



BioMed Central 



